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I investigate energy transfer in a donor-acceptor pair beyond weak system-bath coupling. I identify 
a transition from coherent to incoherent dynamics with increasing temperature, due to multi-phonon 
effects not captured by a standard weak-coupling treatment. The crossover temperature has a 
marked dependence on the degree of spatial correlation between fluctuations experienced at the two 
system sites. For strong correlations, this leads to the possibility of coherence surviving into a high 
temperature regime. 

PACS numbers: 71.35.-y, 03.65.Yz 



Excitation energy transfer is a fundamental process 
common to a wide variety of multi-site (donor-acceptor) 
systems, ranging from those in the solid-state, such as 
crystal impurities 0, 0, [H and quantum dots (QDs) 0, 
IE III i to conjugat ed poly mers and photosynthetic com- 
plexes 0, 0, HO, [Hi H3 ■ In its simplest Forster-Dexter 
(FD) form energy transfer is considered to be incoher- 
ent, resulting from weak donor-acceptor transition-dipolc 
interactions 0. However, recent experimental progress 
in demonstrating quantum coherent energy transfer in a 
number of systems 0, EE EH has highlighted the impor- 
tance of describing transfer dynamics beyond the inco- 
herent regime [13]. Furthermore, such systems are still 
embedded in a larger host matrix, and therefore remain 
susceptible to couplings to their environment fl3j | . The 
resulting interplay between coherent and incoherent pro- 
cesses can fundamentally alter the nature of the transfer 
dynamics, destroying quantum coherent effects and mod- 
ifying the transfer rate. 

To develop a full understanding of any donor-acceptor 
system it is thus crucial to establish the coherent or inco- 
herent nature of the transfer process 

sunn, and to 

ex- 
plore how this changes with variations in donor-acceptor 
separation, system-bath coupling strengths, or tempera- 
ture. For example, the recent demonstration of coherent 
transfer at room temperature in conjugated polymers 
points to the potentially pivotal role played by correlated 
dephasing fluctuations in protecting coherence in these 
systems [?, TfJ 16 1. Furthermore, determining the respec- 
tive roles of coherent and incoherent processes in optimis- 
ing energy transfer efficiency in donor-acceptor networks 
is currently subject to considerable interest 17, ID, 19l ]. 

A number of methods have been put forward to deal 
with the dynamics of coherent energy transfer under the 
influence of an external environment. A popular assump- 
tion is that the system-bath coupling is weak [l9|, l2fj |. 
which leads to Redfield -typ e dynamics involving only 



single-phonon processes [21|. A modified Redfield treat- 
ment, with a broader range of validity, has also been 



suggested [2JA122J. For strong system-bath coupling, FD 
theory has been extended to account for exciton delocal- 
isation over donor and acceptor sites [23|, while the po- 
laron transformation provides a useful tool to investigate 
both the weak and strong coupling regimes 0, 0, . 
The importance of non-Markovian effects has also been 
studied [HH. 

To explore the criteria for coherent energy transfer in 
a donor-acceptor pair in detail, I present here an analyti- 
cal theory of the transfer dynamics capable of interpolat- 
ing between the weak (single-phonon) and strong (multi- 
phonon) system-bath coupling regimes, and correlated to 
independent fluctuations, while still capturing the coher- 
ent dynamics due to the donor-acceptor electronic cou- 
pling. As a main result, I identify a crossover from co- 
herent to incoherent transfer for resonant donor-acceptor 
pairs with increasing temperature, as multi-phonon ef- 
fects become dominant. Such behaviour cannot be de- 
rived from a weak-coupling treatment. I show that the 
critical temperature at which the crossover occurs has 
a pronounced dependence on the degree of correlation 
between fluctuations at each site, leading to the possibil- 
ity of coherent transfer surviving at high temperatures 
in strongly correlated environments, where multi-phonon 
processes are suppressed. 

Consider a pair of two-level systems (j = 1,2) sepa- 
rated by a distance d, with energy transfer interaction 
Vf, coupled linearly to a harmonic environment (ft = 1): 

H = ^ e J \X) ] {X\ + V F (\GX){XG\ + \XG){GX\) 
3=1,2 

k 3=1,2 k 

Here, each system has ground (excited) state \G)j (\X)j) 
and energy Cj, the system-bath couplings are given by g J k , 
and the bath comprises a collection of oscillators of fre- 
quencies Wk and creation (annihilation) operators (6k)- 
Such a model has previously been employed in a range of 



2 



physical settings, see e.g. Refs. @, H, H, 12, 14, 2(J, and mitian conjugate) [l 



could also represent the basic unit of a spin chain [27 1. 
We shall consider system-bath couplings of the form 
9k = l.9k|e lk ' d ^ 2 and flf = |.9k| e_ikd ^ 2 , where position 
dependent phases give rise to correlations between the 
bath-induced fluctuations ex per ienced at each system 
site [3, H, [H . Following Refs. [20|, [Hj , I parameterise the 
transfer interaction as Vf = [Vb / (d / do) 3 } f (d / do) , where 
d = |d|, f(x) = 3^/^/2 crf(a-/v / 2) - 3xe~ x2 / 2 accounts 
for small d corrections to the dipole approximation, and 
do determines when the dipole limit is reached. 



V 2 J2(M"')lPi("),PiW} + h.c), 

I.uj.uj' 



l (1) 

where lu,lu' G {0,±?7}, i] = t 2 + 4V|, A;(w) = 
"fi(uj)/2 + iSi(uj), and we have decomposed the system 
operators as &i(t) = J2 U Pi(^)e~ luJt [Ifl]. The rates 

ll {oj)=e^' 2 ^-{5/^-^/2)5^0)), (2) 



The full Hamiltonian may be decomposed into three 
decoupled subspaccs [\GG),{\XG),\GX)},\XX)}. We 
are interested in the energy transfer dynamics occur- 
ring between the single-excitation states, described by a 
Hamiltonian iJ su b, and we set \XG) = |0), \GX) = |1) to 
identitfy an effective two-state system spanning the sub- 
space [14J . To move into an appropriate basis for the sub- 
sequent perturbation theory, we apply the unitary trans- 
formation H' = e s H snh e- s , where S = |0)(0| J2k( a ^ b i~ 
a k 6 k ) + |l)<l|Ek(A6i " /£&k), with a k = gl/tu k and 
Ac = Sk/^k- As a result, we map our system to the 
polaron-transformed, spin-boson model H' = |<t 2 + 
Vru x + J^k ^k& k fck + Vf {o- x B x + <JyBy), here describing 
the energy transfer dynamics of our donor-acceptor pair 
in the single-excitation subspace, with bath-renormalised 
coupling Vr = BVf [1, 0,12^]. The Pauli matrices, o\ 
(for I = x,y,z), are defined in the basis {|0), |1)}, while 
e = e\ — £2- Bath-induced fluctuations are described by 
B x = (1/2)(B + + B_-2B) smdB y = {-i/2){B_ - B + ), 
where B± = n k -D(±(ak — /3k)) are products of displace- 
ment operators D(±\k) = e ± ( Xkb k~ x * bk ' 24 1 . Assum- 



ing the bath to be in thermal equilibrium at temper- 
ature T, the correlation-dependent renormalisation of 
the coupling strength is determined by B = (B±) = 

, where (3 = 1/fcsT, with 



Boltzmann constant ks- Here, we define a single-site 
spectral density as J(uj) = J2k l.9k| 2 <5(cj — u>k), while 
the factor (1 — F(uj,d)) accounts for the degree of spa- 
tial correlation in the fluctuations at each site. We find 
F(lo, d) = sinc(wcZ/c) in 3-dimensions, assuming k = lo/c, 
and that J(lo) is isotropic. 



are Fourier transforms of the bath correlation functions 

(Bt(T-i(3/2)B x (0)) = (B 2 /2)(e^ (T) +e-^-2),(3) 
(Bl(r-iL3/2)By(0)) = (B 2 /2)(e^ e~^), (4) 

defined in terms of the phonon propagator [3] 



<p(r) 



duj- 



-(1-F(u,d)) 



sinh {(3uj/2) 



(5) 



while Si(uj) = Im J °° dTe lu)T {B\{T)Bi{0)). The most in- 
teresting dynamics of the model can now be explored by 
considering two limiting cases: that of resonant donor 
and acceptor, in which the interplay between coherent 
and incoherent processes is most evident, and that of 
large energy mismatch, often encountered in practise. 

Resonant - The resonant case is of particular impor- 
tance as it demonstrates most clearly how bath-induced 
fluctuations can fundamentally alter the nature of the 
energy transfer process. As we shall see, in the high- 
temperature regime, multi-phonon dephasing effects can 
become dominant, giving rise to a crossover from low- 
temperature coherent dynamics to a high-tempcraturc 
incoherent process. Setting e = 0, we derive from 
Eq. ([27]) a set of Bloch equations governing the time 
evolution of the system state. Taking an initial state 
p'(0) = \XG)(XG\, a single excitation in the donor, 
and transforming out of the polaron frame, we solve 
for the subsequent donor-acceptor population dynamics, 
(o- z )t = Tr s (cr 2/ o(t)), to find 



(<r,)t = e-^X/ 2 



cos ■ 



Mi) . # 

sin — 

£ 2 



(6) 



We now write W = H + H'j, where H'j = 
Vf {o~ x B x + o-yBy) is treated as a perturbation. Pro- 
vided Vr is non-zero, as assumed throughout, this pro- 
cedure is suitable for exploring both single- and multi- 
phonon bath- induced effects on the system dynamics. In 
cases where Vr = 0, we can instead apply a related 
variational approach [29(. Following the standard pro- 
cedure we derive a Markovian master equation de- 
scribing the reduced system dynamics in the polaron- 
transformed Schrodingcr picture (h.c. denotes the Her- 



where £ = y/8V R (2V R + X) - (r x - r 2 ) 2 . Here, 

, (1 + 2JV(2V fl )) 



Ti = 

r 2 = 



V,i 



27x(0) + 7»(2Vk)- 



(1 + N(2V R )) 



2y|7*(0), 



(7) 
(8) 



0U1 



X = 2V*(S v (2Vr) -S v (-2Vr)), and N(w) = (e 
1) . The coherent-incoherent transition thus occurs at 

8Vr(2Vr + A) = (fi - T 2 ) 2 . (9) 
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We shall return to the crossover shortly. First, let's 
consider the dynamics in the weak system-bath cou- 
pling limit. In this case, we expand Eqs. §5§ and (|1]) 
to first order in <p(t), hence keeping only single-phonon 
contributions. We then find T 2 « 0, and thus a 
damping rate fi = ttJ(2V r )(\ - F(2V R ,d))coth/3V R . 
Here, V R = BVp, where we expand B « B [l — 
J °° du>^r-(l — F(u>, d))(coth/3o;/2 — 1)], with vacuum 
term B Q = e" ^ *"^(l-- F ('"><0). From Eq. ©, we find 
that the system performs damped coherent oscillations: 
(a z ) t = e- fl */ 2 [cos (it/2) - (f i/f) sin (ft/2)], with fre- 
quency £ p 



16V* 



rf , where 4Vr > Ti to be consis- 



tent with the original expansion. 

Under what circumstances is a weak-coupling approxi- 
mation appropriate? To address this question it is useful 
to consider an explicit form for the spectral density. As 
an illustration, we choose J(u>) = Alo 3 , describing, for 
example, acoustic pho non induced dephasing with a cou- 
pling strength A (2dl . 0, Hl[. Here, we keep a cutoff 
frequency uj c only in the vacuum terms. From Eq. ((5|) we 
obtain 



<p(T>) 



sech 



tanh (x — t') + tanh (x + t') 



where we scale the time as r' 



2x 

(10) 

7rr//3, and define the 



dimensionless parameters tpo = 2tt 2 A/(3 2 = T 2 /Tq and 
x = nd/c/3 = T/Td- Importantly, we can now identify 
two distinct temperature scales that determine whether 
single-phonon or multi- phonon processes are relevant: 
To, set by A [24J; and Td, which is inversely proportional 
to the separation, and is therefore correlation- dependent. 

Let's consider two cases: (i) when x ^> 1 (T ^> Td, 
weak fluctuation correlations), it can be shown from 
Eq. (jT0|) that ipo alone is suitable as an expansion pa- 
rameter in the bath correlation functions. Hence, ipo <C 1 
defines the single-phonon regime in this case, most easily 
satisfied for large separation d, small A, and low T; (ii) 
when i<1, the strongly-correlated case most easily sat- 
isfied for small d, we expand Eq. (|10[) to second-order in x 
to give (p(f) m <p x 2 ((l - 4tanh 2 r')/3 + tanh 4 t'). Now, 
fox plays the role of an expansion parameter in the cor- 
relation functions, with the single-phonon rate valid for 
ipox 2 <C 1. However, since x is already assumed small 
in this case, it is clear that the single-phonon rate can 
be used at least up to tpo ~ 1, and is therefore valid 
over a much larger range of temperatures and/or cou- 
pling strengths than in case (i). The system is thus far 
better protected from the adverse effects of the environ- 
ment when the fluctuations are highly correlated, and 
hence multi-phonon processes can be suppressed up to 
much higher temperatures. This is shown in the inset 
of Fig. [TJ where the damping rate in Eq. ((6]) is plotted 
against temperature for strong correlations, leading to a 
single-phonon rate valid beyond T/Tq = 1. 



Turning now to the high-tempcrature regime, the rates 
are estimated by expanding Cp(r) about r = 0, where it 
is strongly peaked. Keeping terms up to r 2 order we find 



Ti « 2r 2 



T/2 d2 „2ipa/3„<po(2xcBcYi2x — l)/x 2 

2f3 ^£ e 

7T(fo (x — sech 2 a; tanh x)/x 



(11) 



valid for 2(iV R < 1, with B 2 = e -n 4 x 2 /( Va +n 2 x 2 ) ^ where 
= cj c /7rfcsTo- Further, A w in this limit, hence £ « 
y/16V%- T 2 /4 » iTi/2 in Eq. ©, giving (a z ) t » e" rit . 
Thus, in the high-temperature resonant case, the transfer 
is incoherent, at a rate Ti given in Eq. (jlip . 

The transition between these two regimes, from co- 
herent to incoherent dynamics, is particularly important 
as it allows us to assess up to what critical tempera- 
ture quantum coherent effects might be observed. As 
we have seen, the weak-coupling dynamics is expected to 
be coherent, hence the crossover generally occurs in the 
high-temperature regime, where Eq. (jlip is valid. Then, 
Eq. © simplifies to 8V R = Ti, with the transfer being 
coherent for 8V R > Ti, We use this condition to define a 
critical temperature, T c , above which the dynamics be- 
comes incoherent. From Eq. (jlip we find the implicit 
equation 

T 2 =T -^-^ ; , (12) 



4k B \/Tr{x c - sech x c tanha; c )/x c 

where <p c = T 2 /T$ and x c = T c /T d . It is clear that T c will 
vary in a nontrivial way as a function of donor-acceptor 
separation, through the dependence of Eq. (fT2"|) on x c , 
Bq, and Vp. Again, we consider two limits: (i) as the 
separation becomes large, the "correlation" temperature 
becomes unimportant (Td — > 0) and T c varies only weakly 
with separation through Vp; (ii) at very small separations 
the rates Ti and T 2 tend to zero, while V R — > Vf- Hence, 
in this limit, T c diverges, as we expect; for complete fluc- 
tuation correlation the system behaviour is always coher- 
ent, with no crossover to incoherent dynamics regardless 
of the temperature. 

To illustrate this behaviour, in the main part of Fig.[T] 
we plot the crossover temperature, shown separating the 
coherent and incoherent regimes, as a function of donor- 
acceptor separation. The divergence of T c at small d im- 
plies that coherent dynamics can survive at elevated tem- 
peratures when strong fluctuation correlations suppress 
multi-phonon effects, consistent with recent experimen- 
tal observations 0]. Further, the change in T c behaviour 
from small to large separations can provide information 
on the correlation length of the bath. Specifically, once 
the distance dependence of T c becomes weak there is no 
longer significant correlation between fluctuations at each 
site. 

To give Fig. [1] a relevant experimental context, we 
now estimate Tq and Td for two closely-spaced semi- 
conductor QDs, as realized experimentally in Ref. 0], 
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Coherent 



12 3 4 

d/d =T /T d 

FIG. 1: (Color online) Main: Regimes of resonant energy 
transfer for varying temperature {T /To) and scaled donor- 
acceptor separation. The line T = T c , given by Eq. (|12|l . 
divides the coherent (lower) and incoherent (upper) cases. 
Inset: Resonant damping rate versus T/Tq evaluated numer- 
ically (black, solid line), and by single-phonon (blue, dashed 
line) and high-temperature (red, dotted line) analytical ap- 
proximations. Here, Td/T) = 10. Parameters: uic/ksTo = 5 
and Vb/fesTo = 1. 



included there. In the opposite, high temperature limit 
(ksT ^> e), we again find (<J z )t ~ e _rt , with T = I\ of 
Eq. (dH). 

Summary - I have presented an analytical theory of 
excitation transfer in a correlated environment, show- 
ing that for resonant donor and acceptor, a crossover 
from coherent to incoherent transfer is expected as multi- 
phonon effects begin to dominate. The theory outlined 
here opens up intriguing possibilities for further study of 
the role of coherence in the transfer dynamics of larger 
arrays, such as photosynthetic complexes 17,3,[l3|. For 
example, it enables one to address the important ques- 
tion of how the transfer efficiency changes in such systems 
when crossing from the coherent to incoherent regime. 

I am very grateful to A. Olaya-Castro, S. Bosc, and 
A. M. Stoneham for useful discussions. I am supported 
by the EPSRC, Griffith University, and the Australian 
Research Council Centre for Quantum Computer Tech- 
nology. 



which could be brought into resonance by applying an 
external electric field. Typically, deformation potential 
coupling to acoustic phonons dominates exciton dephas- 
ing in such samples [2(j. A simple model 31 1 allows 
an estimate of A = 0.032 ps 2 in this case, implying 
To ~ 9.6 K. Taking c = 5110 ms" 1 [Io| and a dot sep- 
aration d = 4.5 nm d, we find T d w 2.8 K. Setting 
Vo/ksTo — 1 and d/do = Tq/Tj then implies the reason- 
able values Vq ~ 0.8 meV and do f» 1.3 nm [2(J, respec- 
tively. From Fig. [T] we then obtain a crossover temper- 
ature of T c w 30 K, below which we expect the energy 
transfer dynamics to display signatures of coherence. In 
fact, in Ref. Q temperatures of around 4 — 40 K were ex- 
plored, which should therefore be a promising range over 
which to observe both coherent and incoherent transfer 
dynamics in QD samples. 

Off-resonant - It is also important to examine the dy- 
namics when the donor and acceptor are far off-resonant 
with each other, such that Vp/e -C 1. This can occur 
quite naturally, for example in QD samples due to the 
nature of their growth. Furthermore, the recent weak- 
coupling theory of Ref. [2(J predicts a single transfer rate 
in the off-resonant regime, and thus provides a means 
to assess the validity of our theory in this limit. As 
in the resonant case, we derive a set of Bloch equa- 
tions from Eq. (|27|) . this time expanding the resulting 
expressions to second-order in ~Vp/t. We find system 
dynamics well approximated by {o~ z )t ~ e~ rt — (1 — 
e _rt ) tanh (/3e/2), describing incoherent energy transfer 
from the initially excited donor to the acceptor at a rate 

2 (l+2JV(e)) 



V = V 1 

pling limit of V by retaining only single-phonon terms, 
we find f w (47rV^ 
sistcnt with Ref. 12 C 



(l+w/e)) (^( £ ) + lvi e ))- Taking the weak cou- 



2 ) J(e)(l - F(e, d)) coth Q3e/2), con- 
once renormalisation of Vr has been 
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SUPPLEMENT: DERIVATION OF THE MASTER EQUATION 



In this supplement, I outline the derivation of the master equation (Eq. (1) of the paper) used to describe the donor-acceptor 
energy transfer dynamics in the single-excitation subspace. 

As in the paper, we write the polaron-transformed Hamiltonian within the single excitation subspace {\XG}, \GX}} 

as 

H' = ( -<j z + V R cr x + Wk&iL& k + Vf {<?xB x + cjyBy) , (13) 

k 

where we set \XG) = |0), \GX) = |1), and define a z = |0)(0| - |1)(1|, a x = a + + <r_ = |0)(1| + |1)(0|, and 
<jy — i(<T- — <t+) = (0| — |0) (1|). We now separate the Hamiltonian as H 1 = + H' It where H' Q = H' s + H' B , with 

H' s = ^a z + V a a x , (14) 
H' B = ^ Wk b*6 k , (15) 

k 

while 

H'j = V F {a x B x + a y B y ) , (16) 

will be treated as a perturbation. 

To proceed with deriving the master equation, we first diagonalise the system part of the polaron-transformed 
Hamiltonian, H' s , by applying the rotation U^H'U, where U = e~ l ^ ay ^ 2 and $ = arctan (2Vr/c). This gives 

U f H'U = %j z + ^ w k 6 k 6 k + V F [(singer, + cos &a x )B x + a y B y ] , (17) 



where rj = \J e 2 + 4 V^, cos?9 = e/rj, and sini? = 2Vr/t]. Now, moving into the interaction picture with respect to 
(r]/2)a z + X^k^k^ik^k' wc obtain an interaction Hamiltonian of the form 

H'j(t) = V F (a x (t)B x (t) + a y (t)B y (t)), (18) 

where we decompose the system operators as [l[ 

&i(t)= Yl p i(") e ~ iujt > ( 19 ) 

wS{0,±r;} 



where I = x,y and 



P x (0) = sintfa 2 , 

P x (±77) = COS $(7=p , 

Py(Q) = 0, 

P y (± V ) = ±ia T . (20) 



Furthermore, the bath operators transform as 



B x (t) = i(S + (*) + fl_(i)-2B), (21) 
B y {t) = i(B_(t)-S+(t)), (22) 

where 

B±(t) = n k £>(±(a k - /?k)e^ kt ), (23) 
are written in terms of the (now time-dependent) bath displacement operators D(±Xk) ~ e ziz ( Xkb *~ Xkbk \ 
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We now follow the standard procedure to derive a Markovian master equation [l[ , here governing the dynamics of 
the reduced system density operator p'(t) in the polaron frame. We integrate the von Neumann equation for the joint 
system-bath density operator in the polaron frame interaction picture, x'i{t), then trace over the bath modes. This 
results in an integro-differential equation for the reduced density operator within the interation picture of the form 



(24) 



where we assume factorising initial conditions, x'(0) = p'(0)pb, with pb = e^ H B^ kBT /TrB(e~ HB ^ kBT ) being the 
thermal equilibrium state of the bath, and use TiB(Hj(t)pB) = 0. To perform the Born-Markov approximation, wc 
now make two assumptions. First, that the perturbation of the bath state is weak during the combined system-bath 
evolution, so that we may factorize the joint density operator as Xj(t) = p'j(t)pB at all times. Second, that the 
timescale on which the system evolves appreciably is large compared to the bath memory time tb, allowing us to 
replace p'j(s) by p'j(t) in Eq. ([24]), giving 

=-J*dsTr B ([H>(t), [H'^p'^pb]]) . (25) 

To complete the Markov approximation, we make a change of variable t— s — > r, and take the upper limit of integration 
to infinity, to give 



^ = - T drTrBm^AH'At-rlp'^pB]]) 



(26) 

Substituting Eq. (fl"8|) into Eq. (|26p . transforming out of the interaction picture, and using Eqs. (fl"9]l - (|23|) . leads 
directly to the Markovian master equation given in Eq. (1) of the paper, 

^ = ~W,,At)]-V$ £ (A,( W ')[fi( U ) 1 fl(a)V(t)]+H.c.), (27) 

l,UJ,Uj' 

describing the polaron-transformed Scrodinger picture dynamics on timescales t$ > tb, where tb ~ l/w c at T = 0, 
with uj c being a high-frequency cutoff in the bath spectral density. Here, u>, to' £ {0, ±77}, H.c. denotes the Hermitian 
conjugate, while A;(cj) are one-sided Fourier transforms of the bath correlation functions: 

A ; H = [°° dre™ T (Bj (r)Bj (0)) , (28) 



where (B}(t)Bi(0)) = Tr B (S/ (t)Bi(0) Pb ) are evaluated using Eqs. (2TJ), (E2J), and J23J). Setting Ai(u) = 7i(w)/2 + 
iSi(u>) to separate real and imaginary parts, we see that the rates Ji{oj) may be written precisely as in Eq. (2) of the 
paper, 



= Ai(lu) + Ai(lu)* = / dre <WT (S/(r)Bi(0)), 

J — DC 

/OO 
dre lUT {B\{T - i(3/2)Bi(0)), (29) 
-CSO 

where a change of integration variable r — * r — i/3/2 has been made in the second line 0. 
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